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Abstract. 

Chandra or XMM-Newton observations of quiescent low-mass X-ray binaries can provide important constraints on the 
equation of state of neutron stars. The mass and radius of the neutron star can potentially be determined from fitting a neutron 
star atmosphere model to the observed X-ray spectrum. For a radius measurement it is of critical importance that the distance 
to the source is well constrained since the fractional uncertainty in the radius is at least as large as the fractional uncertainty in 
the distance. Uncertainties in modelling the neutron star atmosphere remain. At this stage it is not yet clear if the soft thermal 
component in the spectra of many quiescent X-ray binaries is variable on timescales too short to be accommodated by the 
cooling neutron star scenario. This can be tested with a long XMM-Newton observation of the neutron star X-ray transient 
Cen X^ in quiescence. With such an observation one can use the Reflection Grating Spectrometer spectrum to constrain the 
interstellar extinction to the source. This removes this parameter from the X-ray spectral fitting of the EPIC pn and MOS 
spectra and allows one to investigate whether the variability observed in the quiescent X-ray spectrum of this source is due to 
variations in the soft thermal spectral component or variations in the power law spectral component coupled with variations 
in Nh- This will test whether the soft thermal component can indeed be due to the hot thermal glow of the neutron star. 
Irrespective of the outcome of such a study, the observed cooling in quiescence in sources for which the crust is significantly 
out of thermal equilibrium with the core due to a prolonged outburst, such as KS 1731-260, seem excellent candidates for 
mass and radius determinations through modelling the observed X-rays with a neutron star atmosphere model (the caveats 
about the source distance and atmosphere modelling do also apply here obviously and presently prevent one from obtaining 
such constraints). Finally, the fact that the soft thermal glow in sources such as SAX J1808.4-3658 and IH 1905-1-000 has not 
been detected in quiescence means that the neutron star cores of these sources must be cold. The most plausible explanation 
seems to be that the neutron stars are more massive than 1.4 M(7, and cool via the direct URCA process. 
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INTRODUCTION 

The majority of low-mass X-ray binaries is found to be 
transient - the so called soft X-ray transients (here, I 
will just refer to them as X-ray transients or just tran- 
sients, see e.g. [1]). The accretion rate in quiescence is 
much reduced or even stopped in these transient systems 
in quiescence. Before the XMM-Newton and Chandra 
satellites were operational only a few neutron star X- 
ray transients had been studied in quiescence (e.g. the 
systems Cen X-A and Aql X-1; [2], [3], [4]). Currently, 
many more systems have been studied in quiescence (see 
e.g. [5], [6], [5], [7], [8], [9]) and Cen X-A and Aql X-1 
have been studied in much more detail than was possible 
before (e.g. [10], [11]). As will become clear later, these 
observations can have a profound impact on an important 
area of astrophysics: determining the neutron star equa- 
tion of state (EoS). This is one of the ultimate goals of 
the study of neutron stars. 

Theoretically, one expects the neutron star to emit 
X-rays even after accretion has stopped ([12]). Due to 



the large heat capacity the temperature of the neutron 
star core is set on timescales of tens of thousands of 
years ([13]). The equilibrium core temperature depends 
on the balance between the heating and the cooling rate 
of the neutron star The heating rate depends on the to- 
tal amount of accreted baryons and the pycnonuclear re- 
actions taking place a few hundred meters deep in the 
crust ([14], [12], [13]). The time-averaged mass accre- 
tion rate in neutron star transients can be derived from 
binary evolution models if the orbital period is known 
(e.g. [15]). The pycnonuclear reactions taking place in 
the neutron star crust are described in [16] and [17]. The 
cooling properties of the neutron star core depend on the 
equation of state (see e.g. [18], [19]). 

The resultant X-ray emission can be modelled by a 
neutron star atmosphere (NSA) model ([20], [21], [22]). 
The NSA models have four free parameters. The neutron 
star distance, mass, radius and temperature. Therefore, in 
theory, an NSA-fit provides means to measure the mass 
and radius of the neutron star and thus constrain the equa- 
tion of state (EoS) of matter at supranuclear densities. In 



practice, spectra of neutron star transients in quiescence 
where the flux is high enough to allow for a spectral 
study, are indeed well-fit by a neutron star atmosphere 
model (NSA). Sometimes, additional emission is present 
at energies above a few keV, often quantified by a power- 
law. When systems are selected for which the distance is 
well known (most notably sources in globular clusters), 
neutron star masses and radii can be determined accu- 
rately (e.g. [23], [24]). Indeed, values for the neutron star 
mass and radius consistent with those of canonical neu- 
tron stars are found, rendering support for the interpreta- 
tion that the soft X-ray spectral component is due to the 
NSA. 

Besides fitting for the neutron star radius and mass in 
the NSA modelling the neutron star mass and possibly 
the neutron star equation of state can be constrained by 
investigating the neutron star temperature. The stringent 
upper limit on the NSA spectral component to the lumi- 
nosity of SAX J1808.4-3658 (< 10%; [7], [25]) and the 
deep hmit on the luminosity and neutron star tempera- 
ture in IH 1905-1-000 ([26], [27]) hint at massive neutron 
stars in these two transients. The reasoning is as follows: 
the upper limits on the thermal spectral component imply 
low core temperatures, which in SAX J1808. 4-3658 and 
IH 1905H-000 are so low that they imply a rapid release 
of the energy produced in the crust. This rapid release of 
energy can only occur via enhanced neutrino emission in 
neutron stars with masses larger than 1.6-1.7 Mq ([18]). 
Neutron stars with masses well above 1.4M0 cannot ex- 
ist for so-called soft equations of state (EoS), in which 
matter at high densities is relatively compressible (e.g., 
due to a meson condensate or a transition between the 
hadron and quark-gluon phases). An important caveat is 
that the heating accretion history of the last several tens 
of thousands of years has to be determined from binary 
evolution models. 

In addition to the short duration transients there exists 
a group of X-ray binaries that return to quiescence after 
a period of activity that lasted several years to more than 
a decade ([5]; for the list of currently known systems 
see [28]). In those sources the neutron star crust is not 
in thermal equiUbrium with the core. The pycnonuclear 
reactions have heated the crust to temperatures higher 
than those related to thermal equilibrium ([6]). 

Here, I briefly discuss these different ways in which 
current and future observations of quiescent low-mass 
X-ray binaries can provide constraints on the neutron 
star equation of state. 

OBSERVED QUIESCENT PROPERTIES 

In Figure 1 I provide an update to the figure from [29]. 
For this I have used results that have been published 
since 2004 and the results from the analysis of an XMM- 



Newton observation of 4U 1608-552 in quiescence (Al- 
tamirano et al. in prep; see Figure 2). As can be seen 
from Figure 1, there seems to be a systematic trend in the 
contribution of the power-law to the total quiescent lu- 
minosity in the 0.5-10 keV band. At luminosities above 
« 10^'^ erg s^' the spectrum can be completely described 
by a power law spectral component. For decreasing lu- 
minosities the power law contribution decreases as well. 
Near luminosities of a few times lO-'^ erg s^' the power 
law contributes w 10%; most of the emission in the 0.5- 
10 keV band can be described by a neutron star atmo- 
sphere or black body model. Also around this luminos- 
ity the correlation between the power-law contribution 
and the luminosity changes to an anti-correlation. Which 
means that at luminosities lower than 10^^ erg s~^ the 
power law becomes progressively more important. In- 
terestingly, this anti-correlation can be explained if the 
power-law contribution remains constant and the soft 
thermal contribution decreases gradually when the lumi- 
nosity changes from 10^^ erg s^' to 10^^ erg s~^. 

RADIUS MEASUREIVIENTS OF 
NEUTRON STARS? 

Heinke et al. [23] and [24] find values for the neutron 
star mass and radius consistent with those of canonical 
neutron stars using neutron star atmosphere model fits 
to X-ray spectra of quiescent X-ray binaries in Glob- 
ular Clusters. However, it is still not certain whether 
the soft thermal spectral component is really due to the 
cooling neutron star. Substantial variations on timescales 
that vary from hundreds of seconds to days and months 
have been observed in the quiescent X-ray luminosity 
of 4U 2129H-47, Aql X-1, IGR J00291-I-5934, flie neu- 
tron star transient in quiescence in NGC 6440 (CXI), 
SAX J1808.4-3658 and Cen X-4 (cf. [30], [11], [31], 
[32], [25], [4]). In the cases where the statistics is best 
i.e. Aql X-1 and Cen X-4, it is clear that the variability 
can not be described by only variations in the power law 
([11], [10]). 

Even though the neutron star mass, radius, distance 
and temperature cannot vary on short timescales, small 
observed changes in the neutron star effective tempera- 
ture can be explained in Ught of an NSA model if an out- 
burst or a type I X-ray burst took place between the ob- 
servations that provide the evidence for variability. Small 
temperature changes can be caused by changes in the 
heat blanketing layer below the neutron star atmosphere 
([33]). The heat blanketing layer consists of ashes of nu- 
clear burning produced in type 1 X-ray bursts and of a 
layer of H and He that remains after an outburst. The 
thickness of the latter layer varies from outburst to out- 
burst. A thicker layer means a higher heat conductivity 
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FIGURE 1. An update of the Figure presented in [29], showing the fractional contribution of the power-law spectral component 
to the 0.5-10 keV quiescent luminosity in the sample of sources bright enough for spectral fitting. The black dots are sources in 
Globular Clusters (from [23]). The data points without an upper error bar denote upper limits. Two points can be noticed. First, the 
quiescent luminosity of some sources is variable. Second, there seems to be a systematic trend in the contribution of the power- 
law to the total quiescent luminosity; at high quiescent luminosities the spectrum is completely described by a power law spectral 
component. Near luminosities of a few times 10"'"' erg s^' the power law contributes < 10%; most of the emission in the 0.5-10 
keV band close to this luminosity can be described by a neutron star atmosphere or black body model. At lower luminosities still 
the power law becomes progressively more important. 
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FIGURE 2. The XMM-Newton pn spectrum of 4U 1608- 
552 in quiescence fitted with a fit-function consisting of a neu- 
tron star atmosphere model plus a power-law component both 
subject to interstellar extinction. To show the contribution of 

the power-law component the normalisation of the power law 
has been set to zero for the plot. A more detailed description 
of the results of our analysis of the XMM-A'evvton data of this 
source will appear in Altamirano et al. (in preparation). 



which implies a higher observed effective temperature 
for a given (unchanged) core temperature. However, the 
luminosity in the soft band in Cen X-4 varies on time 
scales too short to be explained by variations in the thick- 
ness of the H/He layer. 

Currently, the favoured explanation is that the power- 
law spectral component varies in accord with Nh ([34], 
[10]). In this way the thermal spectral component can be 
kept constant. If the soft thermal component varies on 
short time scales there is a problem with the interpreta- 
tion that it arises in the NSA. 

Usage of spectra obtained with the 
XMM-Newton Reflecting Grating 
Spectrometer 

Coupled variations between the power-law index and 
Nh such as found in Aql X-1 by [34] and in Cen 
by [10] can be accommodated in the model that explains 
the power law component as shocked emission from the 
pulsar wind and the infalling material ([35], [34]) if the 
amount of (ionised) matter at the shock region and Nh 
are correlated. The nature of this correlation is unclear. 
Hence, the coupled variability between Nh and the power 
law is rather ad hoc. 

With a long XMM-Newton observation of Cen X-4 



one can determine Nh by measuring the equivalent width 
of the O / K-edge observable in the RGS spectrum (see 
Figure 3; figure from [36]). In this way Nh can be deter- 
mined independently from the broadband spectral fit to 
the EPIC pn and MOS spectra, leaving only the temper- 
ature and normalisation of the soft thermal component 
and the power-law index and normalisation as free pa- 
rameters to explain the variability. 

There are two possible outcomes of such a study. 
First, if it is determined that indeed the soft thermal 
spectral component does not vary and the variability can 
be explained by coupled variations in the power law and 
Nh it would boost the confidence in the masses and radii 
that have been and can be determined from fitting the 
NSA models to the X-ray spectra of quiescent low-mass 
X-ray binaries. Note that unknown uncertainties in the 
neutron star atmosphere models remain. However, in this 
respect it is interesting to note that various atmosphere 
models give similar results (see discussion in [24]). 

The second possible outcome is that the soft thermal 
component varies substantially. This would either mean 
that there exists a currently unidentified mechanism asso- 
ciated with crustal heating that causes the effective tem- 
perature to change on short timescales or that the soft 
X-ray spectral component is caused by another process 
such as residual accretion as considered by [2] and [37]. 

UNDETECTED QUffiSCENT NEUTRON 
STAR X-RAY TRANSIENTS 

If the soft thermal component that is observed in many 
quiescent X-ray transients is not due to the thermal 
glow of the hot neutron star core moderated by the NSA 
this would imply that the cores of these neutron stars 
are so cold that the hot thermal glow is not detectable, 
providing evidence for enhanced coohng mechanisms in 
these neutron stars (cf . [27]for the case of IH 1905h-000). 
We have identified several sources which have not been 
detected in quiescence in X-rays (see Table 1). Note that 
the error on some of the luminosities mentioned in Table 
1 can be of the order of a factor 4 due to uncertainties in 
the distance. 

In the case of IH 1905H-000 the distance is known with 
an accuracy of « 20% as a radius expansion type I X-ray 
burst has been discovered from this source ([45]). The 
extremely deep limit on the quiescent luminosity implies 
that the neutron star core probably has access to the di- 
rect URCA neutrino emission mechanism. This conclu- 
sion can be avoided if neutron star core is not in a steady 
state. For instance if the recurrence time of this transient 
is of the order of ten thousand years (Brown, E.F. pri- 
vate communication). A neutron star transient that has 
been detected in quiescence but for which the soft ther- 
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FIGURE 3. A simulated RGS spectrum for a 250 ks-long XMM-A'ewton observation of Cen X^(from [36]). The black bullets 
with error bars are the rebinned simulated RGS-1 data points, the (blue) open triangles with error bars are the rebinned simulated 
RGS-2 data points, and the (red), drawn line is the model fitted to the RGS-1 data where the O column density is put to zero 
showing the large difference between the data and the model especially between 15-22 A. Most of the narrow lines are due to CCD 
gaps and detector drop outs. Most of the O I K-edge wavelength is not covered by the RGS-2 due to loss of CCD4 in RGS-2. The 
O column and, under the assumption of solar abundances, the Nh can be determined with an accuracy of 8-9% if the broad band 
spectral components are held fixed. Such an observation will allow one to determine if indeed the factor of 2-3 variation in 
necessary to explain the luminosity variations as coupled variations in the power law index and Nh ([10]) is present or not. 



TABLE 1. Neutron star X-ray transients that have not been detected by 
Chandra or XMM-Newton while in quiescence. 



Undetected Sources Luminosity limit (erg s ) Reference 



XMMU J174716.1-281048 


<6.2xl0^2 


[38] 


2S 1803-245 


<2.8xl0^2 


[39] 


EXO 1747-214 


<7xl0^' 


[40] 


Terzan 1 


<lxlO^^ * 


[41] 


XTEJ175 1-305 


<0.2-2xl0^2 


[42] 


XTE J 1807-294 


<4xl03^ 


[43] 


IH 1905+000 


<2xl03° 


[26], [27] 



* but see [44] 



mal spectral component remains so far undetectable is 
SAX J1808. 4-3658. A detailed discussion of that source 
is given by Heinke et al. in this volume. 

LONG-DURATION NEUTRON STAR 
X-RAY TRANSffiNTS IN QUIESCENCE 

There exists a group of transients for which the out- 
burst(s) last so long that after an outburst ends the crust 
is no longer in thermal equilibrium with the core. When 
the accretion rate decreases/stops the cooling neutron 
star crust can be observed (e.g. [5], [6], [46]). The two 
sources that have been observed in detail during this 
cooling phase are KS 1731-260 and MXB 1658-29. The 
latest results on these sources have been presented by 
[46]. Note that even though these sources are also plot- 
ted in Figure 1, their spectra can be completely described 
by an NSA: only upper limits on a power law contribu- 
tion have been derived. Furthermore, even if the soft ther- 
mal component in Cen turns out to be variable this 
does not necessarily imply that the thermal component 
observed in these long-duration transients is of the same 
origin as that in the other (short-duration) transients. 
Therefore, presently, NSA fits to the spectra of such 
long-duration transients just after an outburst seem an 
excellent way to determining the neutron star mass and 
radius. Unfortunately, KS 1731-260 and MXB 1658-29 
suffer from a rather high Nh and their distance is large 
Umiting the source flux in the soft X-ray bands. Fur- 
thermore, even though their distances have been derived 
via radius expansion bursts, the distance uncertainty is 
still too large to provide stringent constraints on the neu- 
tron star radius in these two sources. A neutron star low- 
mass X-ray binary in a nearby globular cluster that has 
been accreting persistently for several years turning to 
quiescence would provide an ideal source to follow with 
Chandra and or XMM-Newton. 

Cackett et al. [46] present evidence for direct URCA 
neutrino emission processes taking place in the cores of 
KS 1731-260 and MXB 1658-29 and for a high termal 
conductivity of their crusts. Recent results of [47], how- 
ever, show that the observations of KS 1731-260 can be 
explained with modified URCA and a "normal" crustal 
thermal conductivity (MXB 1658-29 has not been mod- 
elled in detail by these authors yet). 

CONCLUSIONS 

There are several potential ways in which (X-ray) obser- 
vations of quiescent low-mass X-ray binaries can pro- 
vide constraints on the neutron star equation of state via 
mass and radius determinations, but also via the derived 



cooling properties. Chandra and XMM-A'ewton obser- 
vations of a future nearby transient that had an outburst 
long enough to heat the crust significantly out of ther- 
mal equilibrium with the core should provide constraints 
on the neutron star mass and radius if the source dis- 
tance can be determined accurately. For the more com- 
mon transients, i.e. those who returned to quiescence af- 
ter a weeks to months-long outburst, the observed spec- 
tral variability should be investigate further before firm 
conclusions about the neutron star radius and mass can 
be drawn from fitting neutron star atmosphere models to 
the quiescent X-ray spectra. A long XMM-Newton ob- 
servation of e.g. Cen X^ would allow the RGS spec- 
trum to be extracted and to investigate if indeed the soft 
thermal spectral component is stationary. If so, and if 
the distance to the system is accurately known the neu- 
tron star mass and radius can be determined from a NSA 
model fit to the quiescent spectra. Finally, the low qui- 
escent (thermal) luminosities of several sources (most 
notably IH 1905H-000 and SAX J1808.4-3658) provide 
evidence for the presence of neutrino emissivities larger 
than can be provided by the modified URCA process. 
Presently, it is thought that only neutron stars more mas- 
sive than the canonical 1.4 M© might have access to 
these direct URCA processes. 
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